\(\int \frac {1}{(3+3 \sin (e+f x))^{5/2} (c+d \sin (e+f x))^{3/2}} \, dx\) [604]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (warning: unable to verify)
   Maple [B] (verified)
   Fricas [B] (verification not implemented)
   Sympy [F(-1)]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 29, antiderivative size = 256 \[ \int \frac {1}{(3+3 \sin (e+f x))^{5/2} (c+d \sin (e+f x))^{3/2}} \, dx=-\frac {\left (c^2-6 c d+25 d^2\right ) \text {arctanh}\left (\frac {\sqrt {\frac {3}{2}} \sqrt {c-d} \cos (e+f x)}{\sqrt {3+3 \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}\right )}{48 \sqrt {6} (c-d)^{7/2} f}-\frac {\cos (e+f x)}{4 (c-d) f (3+3 \sin (e+f x))^{5/2} \sqrt {c+d \sin (e+f x)}}-\frac {(3 c-13 d) \cos (e+f x)}{48 (c-d)^2 f (3+3 \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}}-\frac {(c-7 d) d (3 c+7 d) \cos (e+f x)}{144 (c-d)^3 (c+d) f \sqrt {3+3 \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \]

[Out]

-3/32*(c^2-6*c*d+25*d^2)*arctanh(1/2*cos(f*x+e)*a^(1/2)*(c-d)^(1/2)*2^(1/2)/(a+a*sin(f*x+e))^(1/2)/(c+d*sin(f*
x+e))^(1/2))/a^(5/2)/(c-d)^(7/2)/f*2^(1/2)-1/4*cos(f*x+e)/(c-d)/f/(a+a*sin(f*x+e))^(5/2)/(c+d*sin(f*x+e))^(1/2
)-1/16*(3*c-13*d)*cos(f*x+e)/a/(c-d)^2/f/(a+a*sin(f*x+e))^(3/2)/(c+d*sin(f*x+e))^(1/2)-1/16*(c-7*d)*d*(3*c+7*d
)*cos(f*x+e)/a^2/(c-d)^3/(c+d)/f/(a+a*sin(f*x+e))^(1/2)/(c+d*sin(f*x+e))^(1/2)

Rubi [A] (verified)

Time = 0.58 (sec) , antiderivative size = 270, normalized size of antiderivative = 1.05, number of steps used = 6, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.207, Rules used = {2845, 3057, 3063, 12, 2861, 214} \[ \int \frac {1}{(3+3 \sin (e+f x))^{5/2} (c+d \sin (e+f x))^{3/2}} \, dx=-\frac {3 \left (c^2-6 c d+25 d^2\right ) \text {arctanh}\left (\frac {\sqrt {a} \sqrt {c-d} \cos (e+f x)}{\sqrt {2} \sqrt {a \sin (e+f x)+a} \sqrt {c+d \sin (e+f x)}}\right )}{16 \sqrt {2} a^{5/2} f (c-d)^{7/2}}-\frac {d (c-7 d) (3 c+7 d) \cos (e+f x)}{16 a^2 f (c-d)^3 (c+d) \sqrt {a \sin (e+f x)+a} \sqrt {c+d \sin (e+f x)}}-\frac {(3 c-13 d) \cos (e+f x)}{16 a f (c-d)^2 (a \sin (e+f x)+a)^{3/2} \sqrt {c+d \sin (e+f x)}}-\frac {\cos (e+f x)}{4 f (c-d) (a \sin (e+f x)+a)^{5/2} \sqrt {c+d \sin (e+f x)}} \]

[In]

Int[1/((a + a*Sin[e + f*x])^(5/2)*(c + d*Sin[e + f*x])^(3/2)),x]

[Out]

(-3*(c^2 - 6*c*d + 25*d^2)*ArcTanh[(Sqrt[a]*Sqrt[c - d]*Cos[e + f*x])/(Sqrt[2]*Sqrt[a + a*Sin[e + f*x]]*Sqrt[c
 + d*Sin[e + f*x]])])/(16*Sqrt[2]*a^(5/2)*(c - d)^(7/2)*f) - Cos[e + f*x]/(4*(c - d)*f*(a + a*Sin[e + f*x])^(5
/2)*Sqrt[c + d*Sin[e + f*x]]) - ((3*c - 13*d)*Cos[e + f*x])/(16*a*(c - d)^2*f*(a + a*Sin[e + f*x])^(3/2)*Sqrt[
c + d*Sin[e + f*x]]) - ((c - 7*d)*d*(3*c + 7*d)*Cos[e + f*x])/(16*a^2*(c - d)^3*(c + d)*f*Sqrt[a + a*Sin[e + f
*x]]*Sqrt[c + d*Sin[e + f*x]])

Rule 12

Int[(a_)*(u_), x_Symbol] :> Dist[a, Int[u, x], x] /; FreeQ[a, x] &&  !MatchQ[u, (b_)*(v_) /; FreeQ[b, x]]

Rule 214

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x/Rt[-a/b, 2]], x] /; FreeQ[{a, b},
x] && NegQ[a/b]

Rule 2845

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Sim
p[b^2*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x])^(n + 1)/(a*f*(2*m + 1)*(b*c - a*d))), x] + Dis
t[1/(a*(2*m + 1)*(b*c - a*d)), Int[(a + b*Sin[e + f*x])^(m + 1)*(c + d*Sin[e + f*x])^n*Simp[b*c*(m + 1) - a*d*
(2*m + n + 2) + b*d*(m + n + 2)*Sin[e + f*x], x], x], x] /; FreeQ[{a, b, c, d, e, f, n}, x] && NeQ[b*c - a*d,
0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0] && LtQ[m, -1] &&  !GtQ[n, 0] && (IntegersQ[2*m, 2*n] || (IntegerQ
[m] && EqQ[c, 0]))

Rule 2861

Int[1/(Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]*Sqrt[(c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]]), x_Symbol] :> D
ist[-2*(a/f), Subst[Int[1/(2*b^2 - (a*c - b*d)*x^2), x], x, b*(Cos[e + f*x]/(Sqrt[a + b*Sin[e + f*x]]*Sqrt[c +
 d*Sin[e + f*x]]))], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 -
 d^2, 0]

Rule 3057

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])*((c_.) + (d_.)*sin[(e_
.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp[b*(A*b - a*B)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*
x])^(n + 1)/(a*f*(2*m + 1)*(b*c - a*d))), x] + Dist[1/(a*(2*m + 1)*(b*c - a*d)), Int[(a + b*Sin[e + f*x])^(m +
 1)*(c + d*Sin[e + f*x])^n*Simp[B*(a*c*m + b*d*(n + 1)) + A*(b*c*(m + 1) - a*d*(2*m + n + 2)) + d*(A*b - a*B)*
(m + n + 2)*Sin[e + f*x], x], x], x] /; FreeQ[{a, b, c, d, e, f, A, B, n}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2
- b^2, 0] && NeQ[c^2 - d^2, 0] && LtQ[m, -2^(-1)] &&  !GtQ[n, 0] && IntegerQ[2*m] && (IntegerQ[2*n] || EqQ[c,
0])

Rule 3063

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])*((c_.) + (d_.)*sin[(e_
.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp[(B*c - A*d)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x]
)^(n + 1)/(f*(n + 1)*(c^2 - d^2))), x] + Dist[1/(b*(n + 1)*(c^2 - d^2)), Int[(a + b*Sin[e + f*x])^m*(c + d*Sin
[e + f*x])^(n + 1)*Simp[A*(a*d*m + b*c*(n + 1)) - B*(a*c*m + b*d*(n + 1)) + b*(B*c - A*d)*(m + n + 2)*Sin[e +
f*x], x], x], x] /; FreeQ[{a, b, c, d, e, f, A, B, m}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2
 - d^2, 0] && LtQ[n, -1] && (IntegerQ[n] || EqQ[m + 1/2, 0])

Rubi steps \begin{align*} \text {integral}& = -\frac {\cos (e+f x)}{4 (c-d) f (a+a \sin (e+f x))^{5/2} \sqrt {c+d \sin (e+f x)}}-\frac {\int \frac {-\frac {3}{2} a (c-3 d)-2 a d \sin (e+f x)}{(a+a \sin (e+f x))^{3/2} (c+d \sin (e+f x))^{3/2}} \, dx}{4 a^2 (c-d)} \\ & = -\frac {\cos (e+f x)}{4 (c-d) f (a+a \sin (e+f x))^{5/2} \sqrt {c+d \sin (e+f x)}}-\frac {(3 c-13 d) \cos (e+f x)}{16 a (c-d)^2 f (a+a \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}}+\frac {\int \frac {\frac {1}{4} a^2 \left (3 c^2-12 c d+49 d^2\right )+\frac {1}{2} a^2 (3 c-13 d) d \sin (e+f x)}{\sqrt {a+a \sin (e+f x)} (c+d \sin (e+f x))^{3/2}} \, dx}{8 a^4 (c-d)^2} \\ & = -\frac {\cos (e+f x)}{4 (c-d) f (a+a \sin (e+f x))^{5/2} \sqrt {c+d \sin (e+f x)}}-\frac {(3 c-13 d) \cos (e+f x)}{16 a (c-d)^2 f (a+a \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}}-\frac {(c-7 d) d (3 c+7 d) \cos (e+f x)}{16 a^2 (c-d)^3 (c+d) f \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}-\frac {\int -\frac {3 a^3 (c+d) \left (c^2-6 c d+25 d^2\right )}{8 \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx}{4 a^5 (c-d)^3 (c+d)} \\ & = -\frac {\cos (e+f x)}{4 (c-d) f (a+a \sin (e+f x))^{5/2} \sqrt {c+d \sin (e+f x)}}-\frac {(3 c-13 d) \cos (e+f x)}{16 a (c-d)^2 f (a+a \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}}-\frac {(c-7 d) d (3 c+7 d) \cos (e+f x)}{16 a^2 (c-d)^3 (c+d) f \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}+\frac {\left (3 \left (c^2-6 c d+25 d^2\right )\right ) \int \frac {1}{\sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx}{32 a^2 (c-d)^3} \\ & = -\frac {\cos (e+f x)}{4 (c-d) f (a+a \sin (e+f x))^{5/2} \sqrt {c+d \sin (e+f x)}}-\frac {(3 c-13 d) \cos (e+f x)}{16 a (c-d)^2 f (a+a \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}}-\frac {(c-7 d) d (3 c+7 d) \cos (e+f x)}{16 a^2 (c-d)^3 (c+d) f \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}-\frac {\left (3 \left (c^2-6 c d+25 d^2\right )\right ) \text {Subst}\left (\int \frac {1}{2 a^2-(a c-a d) x^2} \, dx,x,\frac {a \cos (e+f x)}{\sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}\right )}{16 a (c-d)^3 f} \\ & = -\frac {3 \left (c^2-6 c d+25 d^2\right ) \text {arctanh}\left (\frac {\sqrt {a} \sqrt {c-d} \cos (e+f x)}{\sqrt {2} \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}\right )}{16 \sqrt {2} a^{5/2} (c-d)^{7/2} f}-\frac {\cos (e+f x)}{4 (c-d) f (a+a \sin (e+f x))^{5/2} \sqrt {c+d \sin (e+f x)}}-\frac {(3 c-13 d) \cos (e+f x)}{16 a (c-d)^2 f (a+a \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}}-\frac {(c-7 d) d (3 c+7 d) \cos (e+f x)}{16 a^2 (c-d)^3 (c+d) f \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \\ \end{align*}

Mathematica [A] (warning: unable to verify)

Time = 9.20 (sec) , antiderivative size = 465, normalized size of antiderivative = 1.82 \[ \int \frac {1}{(3+3 \sin (e+f x))^{5/2} (c+d \sin (e+f x))^{3/2}} \, dx=\frac {\left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^4 \left (\frac {\left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right ) \left (-14 c^3+25 c^2 d+56 c d^2+113 d^3+d \left (3 c^2-14 c d-49 d^2\right ) \cos (2 (e+f x))+\left (-6 c^3+14 c^2 d+62 c d^2+170 d^3\right ) \sin (e+f x)\right )}{(c+d) \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^3}+\frac {3 \left (c^2-6 c d+25 d^2\right ) \left (\log \left (1+\tan \left (\frac {1}{2} (e+f x)\right )\right )-\log \left (c-d+2 \sqrt {c-d} \sqrt {\frac {1}{1+\cos (e+f x)}} \sqrt {c+d \sin (e+f x)}+(-c+d) \tan \left (\frac {1}{2} (e+f x)\right )\right )\right )}{\frac {\sec ^2\left (\frac {1}{2} (e+f x)\right )}{2+2 \tan \left (\frac {1}{2} (e+f x)\right )}-\frac {-\frac {1}{2} (c-d) \sec ^2\left (\frac {1}{2} (e+f x)\right )+\frac {\sqrt {c-d} \left (\frac {1}{1+\cos (e+f x)}\right )^{3/2} (d+d \cos (e+f x)+c \sin (e+f x))}{\sqrt {c+d \sin (e+f x)}}}{c-d+2 \sqrt {c-d} \sqrt {\frac {1}{1+\cos (e+f x)}} \sqrt {c+d \sin (e+f x)}+(-c+d) \tan \left (\frac {1}{2} (e+f x)\right )}}\right )}{288 \sqrt {3} (c-d)^3 f (1+\sin (e+f x))^{5/2} \sqrt {c+d \sin (e+f x)}} \]

[In]

Integrate[1/((3 + 3*Sin[e + f*x])^(5/2)*(c + d*Sin[e + f*x])^(3/2)),x]

[Out]

((Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^4*(((Cos[(e + f*x)/2] - Sin[(e + f*x)/2])*(-14*c^3 + 25*c^2*d + 56*c*d^
2 + 113*d^3 + d*(3*c^2 - 14*c*d - 49*d^2)*Cos[2*(e + f*x)] + (-6*c^3 + 14*c^2*d + 62*c*d^2 + 170*d^3)*Sin[e +
f*x]))/((c + d)*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^3) + (3*(c^2 - 6*c*d + 25*d^2)*(Log[1 + Tan[(e + f*x)/2]
] - Log[c - d + 2*Sqrt[c - d]*Sqrt[(1 + Cos[e + f*x])^(-1)]*Sqrt[c + d*Sin[e + f*x]] + (-c + d)*Tan[(e + f*x)/
2]]))/(Sec[(e + f*x)/2]^2/(2 + 2*Tan[(e + f*x)/2]) - (-1/2*((c - d)*Sec[(e + f*x)/2]^2) + (Sqrt[c - d]*((1 + C
os[e + f*x])^(-1))^(3/2)*(d + d*Cos[e + f*x] + c*Sin[e + f*x]))/Sqrt[c + d*Sin[e + f*x]])/(c - d + 2*Sqrt[c -
d]*Sqrt[(1 + Cos[e + f*x])^(-1)]*Sqrt[c + d*Sin[e + f*x]] + (-c + d)*Tan[(e + f*x)/2]))))/(288*Sqrt[3]*(c - d)
^3*f*(1 + Sin[e + f*x])^(5/2)*Sqrt[c + d*Sin[e + f*x]])

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(4112\) vs. \(2(235)=470\).

Time = 5.11 (sec) , antiderivative size = 4113, normalized size of antiderivative = 16.07

method result size
default \(\text {Expression too large to display}\) \(4113\)

[In]

int(1/(a+a*sin(f*x+e))^(5/2)/(c+d*sin(f*x+e))^(3/2),x,method=_RETURNVERBOSE)

[Out]

1/16/f*(-85*cos(f*x+e)*sin(f*x+e)*(2*c-2*d)^(1/2)*d^3-11*cos(f*x+e)*(2*c-2*d)^(1/2)*c^2*d-35*cos(f*x+e)*(2*c-2
*d)^(1/2)*c*d^2+6*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+
e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+s
in(f*x+e)))*c^3+150*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*
x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1
+sin(f*x+e)))*d^3+6*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*
x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1
+sin(f*x+e)))*sin(f*x+e)*cos(f*x+e)*c^3+150*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1
/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*
x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*sin(f*x+e)*cos(f*x+e)*d^3-30*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x
+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+
sin(f*x+e)))*c^2*d*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+114*ln(2*((2*c-2*d)^(1/2)*2^(1/2
)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/
(-cos(f*x+e)+1+sin(f*x+e)))*c*d^2*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)-30*2^(1/2)*((c+d*
sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*
x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*c^2*cos(f*x+e)*d+114
*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1
))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*c*cos
(f*x+e)*d^2+15*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))
/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(
f*x+e)))*c^2*cos(f*x+e)^2*d-57*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*(
(c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-c
os(f*x+e)+1+sin(f*x+e)))*c*cos(f*x+e)^2*d^2+15*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)
^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos
(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*c^2*cos(f*x+e)^3*d-57*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)
*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*
cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*c*cos(f*x+e)^3*d^2+7*cos(f*x+e)*c^3*(2*c-2*d)^(1/2)-4
*sin(f*x+e)*c^3*(2*c-2*d)^(1/2)-4*sin(f*x+e)*d^3*(2*c-2*d)^(1/2)-75*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^
(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+
e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*cos(f*x+e)^3*d^3-3*2^(1/2)*((c+d*sin(f*x+e))/(co
s(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+
e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*c^3*cos(f*x+e)^3+150*2^(1/2)*((c+d*
sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*
x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*cos(f*x+e)*d^3-3*2^(
1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(
1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*c^3*cos(f
*x+e)^2-75*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(co
s(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+
e)))*cos(f*x+e)^2*d^3+6*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*si
n(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+
e)+1+sin(f*x+e)))*c^3*cos(f*x+e)+6*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f
*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*c^3*2^(1/2)*((c+d*s
in(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+150*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1)
)^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*d^3*2^
(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)-30*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln
(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos
(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*c^2*d+114*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2
)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e)+c
*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*c*d^2-30*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(
1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e)+c*sin(f*x+e)-d*sin(f*x+e
)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*sin(f*x+e)*cos(f*x+e)*c^2*d+114*2^(1/2)*((c+d*sin
(f*x+e))/(cos(f*x+e)+1))^(1/2)*ln(2*((2*c-2*d)^(1/2)*2^(1/2)*((c+d*sin(f*x+e))/(cos(f*x+e)+1))^(1/2)*sin(f*x+e
)+c*sin(f*x+e)-d*sin(f*x+e)+c*cos(f*x+e)-d*cos(f*x+e)-c+d)/(-cos(f*x+e)+1+sin(f*x+e)))*sin(f*x+e)*cos(f*x+e)*c
*d^2-4*(2*c-2*d)^(1/2)*c^2*d-4*(2*c-2*d)^(1/2)*c*d^2-31*cos(f*x+e)*sin(f*x+e)*(2*c-2*d)^(1/2)*c*d^2-10*cos(f*x
+e)^2*c^2*d*(2*c-2*d)^(1/2)-17*cos(f*x+e)^2*c*d^2*(2*c-2*d)^(1/2)+3*cos(f*x+e)*sin(f*x+e)*c^3*(2*c-2*d)^(1/2)+
4*sin(f*x+e)*c^2*d*(2*c-2*d)^(1/2)+4*sin(f*x+e)*c*d^2*(2*c-2*d)^(1/2)+49*(2*c-2*d)^(1/2)*cos(f*x+e)^3*d^3-36*(
2*c-2*d)^(1/2)*cos(f*x+e)^2*d^3+4*c^3*(2*c-2*d)^(1/2)+3*cos(f*x+e)^2*sin(f*x+e)*c^2*d*(2*c-2*d)^(1/2)-14*cos(f
*x+e)^2*sin(f*x+e)*c*d^2*(2*c-2*d)^(1/2)-3*cos(f*x+e)^3*c^2*d*(2*c-2*d)^(1/2)+14*cos(f*x+e)^3*c*d^2*(2*c-2*d)^
(1/2)-49*cos(f*x+e)^2*sin(f*x+e)*d^3*(2*c-2*d)^(1/2)+4*(2*c-2*d)^(1/2)*d^3+3*cos(f*x+e)^2*c^3*(2*c-2*d)^(1/2)-
81*cos(f*x+e)*(2*c-2*d)^(1/2)*d^3-7*cos(f*x+e)*sin(f*x+e)*c^2*d*(2*c-2*d)^(1/2))*(c+d*sin(f*x+e))^(1/2)/(d*cos
(f*x+e)^3+cos(f*x+e)^2*d*sin(f*x+e)+c*cos(f*x+e)^2-cos(f*x+e)*sin(f*x+e)*c+2*cos(f*x+e)^2*d-d*sin(f*x+e)*cos(f
*x+e)-c*cos(f*x+e)-2*c*sin(f*x+e)-d*cos(f*x+e)-2*d*sin(f*x+e)-2*c-2*d)/(a*(sin(f*x+e)+1))^(1/2)/a^2/(c+d)/(2*c
-2*d)^(1/2)/(c-d)^3

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1375 vs. \(2 (235) = 470\).

Time = 0.97 (sec) , antiderivative size = 2984, normalized size of antiderivative = 11.66 \[ \int \frac {1}{(3+3 \sin (e+f x))^{5/2} (c+d \sin (e+f x))^{3/2}} \, dx=\text {Too large to display} \]

[In]

integrate(1/(a+a*sin(f*x+e))^(5/2)/(c+d*sin(f*x+e))^(3/2),x, algorithm="fricas")

[Out]

[-1/128*(3*((c^3*d - 5*c^2*d^2 + 19*c*d^3 + 25*d^4)*cos(f*x + e)^4 + 4*c^4 - 16*c^3*d + 56*c^2*d^2 + 176*c*d^3
 + 100*d^4 - (c^4 - 3*c^3*d + 9*c^2*d^2 + 63*c*d^3 + 50*d^4)*cos(f*x + e)^3 - (3*c^4 - 10*c^3*d + 32*c^2*d^2 +
 170*c*d^3 + 125*d^4)*cos(f*x + e)^2 + 2*(c^4 - 4*c^3*d + 14*c^2*d^2 + 44*c*d^3 + 25*d^4)*cos(f*x + e) + (4*c^
4 - 16*c^3*d + 56*c^2*d^2 + 176*c*d^3 + 100*d^4 - (c^3*d - 5*c^2*d^2 + 19*c*d^3 + 25*d^4)*cos(f*x + e)^3 - (c^
4 - 2*c^3*d + 4*c^2*d^2 + 82*c*d^3 + 75*d^4)*cos(f*x + e)^2 + 2*(c^4 - 4*c^3*d + 14*c^2*d^2 + 44*c*d^3 + 25*d^
4)*cos(f*x + e))*sin(f*x + e))*sqrt(2*a*c - 2*a*d)*log(((a*c^2 - 14*a*c*d + 17*a*d^2)*cos(f*x + e)^3 - 4*a*c^2
 - 8*a*c*d - 4*a*d^2 - (13*a*c^2 - 22*a*c*d - 3*a*d^2)*cos(f*x + e)^2 + 4*((c - 3*d)*cos(f*x + e)^2 - (3*c - d
)*cos(f*x + e) + ((c - 3*d)*cos(f*x + e) + 4*c - 4*d)*sin(f*x + e) - 4*c + 4*d)*sqrt(2*a*c - 2*a*d)*sqrt(a*sin
(f*x + e) + a)*sqrt(d*sin(f*x + e) + c) - 2*(9*a*c^2 - 14*a*c*d + 9*a*d^2)*cos(f*x + e) - (4*a*c^2 + 8*a*c*d +
 4*a*d^2 - (a*c^2 - 14*a*c*d + 17*a*d^2)*cos(f*x + e)^2 - 2*(7*a*c^2 - 18*a*c*d + 7*a*d^2)*cos(f*x + e))*sin(f
*x + e))/(cos(f*x + e)^3 + 3*cos(f*x + e)^2 + (cos(f*x + e)^2 - 2*cos(f*x + e) - 4)*sin(f*x + e) - 2*cos(f*x +
 e) - 4)) + 8*(4*c^4 - 8*c^3*d + 8*c*d^3 - 4*d^4 - (3*c^3*d - 17*c^2*d^2 - 35*c*d^3 + 49*d^4)*cos(f*x + e)^3 +
 (3*c^4 - 13*c^3*d - 7*c^2*d^2 - 19*c*d^3 + 36*d^4)*cos(f*x + e)^2 + (7*c^4 - 18*c^3*d - 24*c^2*d^2 - 46*c*d^3
 + 81*d^4)*cos(f*x + e) - (4*c^4 - 8*c^3*d + 8*c*d^3 - 4*d^4 - (3*c^3*d - 17*c^2*d^2 - 35*c*d^3 + 49*d^4)*cos(
f*x + e)^2 - (3*c^4 - 10*c^3*d - 24*c^2*d^2 - 54*c*d^3 + 85*d^4)*cos(f*x + e))*sin(f*x + e))*sqrt(a*sin(f*x +
e) + a)*sqrt(d*sin(f*x + e) + c))/((a^3*c^5*d - 3*a^3*c^4*d^2 + 2*a^3*c^3*d^3 + 2*a^3*c^2*d^4 - 3*a^3*c*d^5 +
a^3*d^6)*f*cos(f*x + e)^4 - (a^3*c^6 - a^3*c^5*d - 4*a^3*c^4*d^2 + 6*a^3*c^3*d^3 + a^3*c^2*d^4 - 5*a^3*c*d^5 +
 2*a^3*d^6)*f*cos(f*x + e)^3 - (3*a^3*c^6 - 4*a^3*c^5*d - 9*a^3*c^4*d^2 + 16*a^3*c^3*d^3 + a^3*c^2*d^4 - 12*a^
3*c*d^5 + 5*a^3*d^6)*f*cos(f*x + e)^2 + 2*(a^3*c^6 - 2*a^3*c^5*d - a^3*c^4*d^2 + 4*a^3*c^3*d^3 - a^3*c^2*d^4 -
 2*a^3*c*d^5 + a^3*d^6)*f*cos(f*x + e) + 4*(a^3*c^6 - 2*a^3*c^5*d - a^3*c^4*d^2 + 4*a^3*c^3*d^3 - a^3*c^2*d^4
- 2*a^3*c*d^5 + a^3*d^6)*f - ((a^3*c^5*d - 3*a^3*c^4*d^2 + 2*a^3*c^3*d^3 + 2*a^3*c^2*d^4 - 3*a^3*c*d^5 + a^3*d
^6)*f*cos(f*x + e)^3 + (a^3*c^6 - 7*a^3*c^4*d^2 + 8*a^3*c^3*d^3 + 3*a^3*c^2*d^4 - 8*a^3*c*d^5 + 3*a^3*d^6)*f*c
os(f*x + e)^2 - 2*(a^3*c^6 - 2*a^3*c^5*d - a^3*c^4*d^2 + 4*a^3*c^3*d^3 - a^3*c^2*d^4 - 2*a^3*c*d^5 + a^3*d^6)*
f*cos(f*x + e) - 4*(a^3*c^6 - 2*a^3*c^5*d - a^3*c^4*d^2 + 4*a^3*c^3*d^3 - a^3*c^2*d^4 - 2*a^3*c*d^5 + a^3*d^6)
*f)*sin(f*x + e)), -1/64*(3*((c^3*d - 5*c^2*d^2 + 19*c*d^3 + 25*d^4)*cos(f*x + e)^4 + 4*c^4 - 16*c^3*d + 56*c^
2*d^2 + 176*c*d^3 + 100*d^4 - (c^4 - 3*c^3*d + 9*c^2*d^2 + 63*c*d^3 + 50*d^4)*cos(f*x + e)^3 - (3*c^4 - 10*c^3
*d + 32*c^2*d^2 + 170*c*d^3 + 125*d^4)*cos(f*x + e)^2 + 2*(c^4 - 4*c^3*d + 14*c^2*d^2 + 44*c*d^3 + 25*d^4)*cos
(f*x + e) + (4*c^4 - 16*c^3*d + 56*c^2*d^2 + 176*c*d^3 + 100*d^4 - (c^3*d - 5*c^2*d^2 + 19*c*d^3 + 25*d^4)*cos
(f*x + e)^3 - (c^4 - 2*c^3*d + 4*c^2*d^2 + 82*c*d^3 + 75*d^4)*cos(f*x + e)^2 + 2*(c^4 - 4*c^3*d + 14*c^2*d^2 +
 44*c*d^3 + 25*d^4)*cos(f*x + e))*sin(f*x + e))*sqrt(-2*a*c + 2*a*d)*arctan(1/4*sqrt(-2*a*c + 2*a*d)*sqrt(a*si
n(f*x + e) + a)*((c - 3*d)*sin(f*x + e) - 3*c + d)*sqrt(d*sin(f*x + e) + c)/((a*c*d - a*d^2)*cos(f*x + e)*sin(
f*x + e) + (a*c^2 - a*c*d)*cos(f*x + e))) + 4*(4*c^4 - 8*c^3*d + 8*c*d^3 - 4*d^4 - (3*c^3*d - 17*c^2*d^2 - 35*
c*d^3 + 49*d^4)*cos(f*x + e)^3 + (3*c^4 - 13*c^3*d - 7*c^2*d^2 - 19*c*d^3 + 36*d^4)*cos(f*x + e)^2 + (7*c^4 -
18*c^3*d - 24*c^2*d^2 - 46*c*d^3 + 81*d^4)*cos(f*x + e) - (4*c^4 - 8*c^3*d + 8*c*d^3 - 4*d^4 - (3*c^3*d - 17*c
^2*d^2 - 35*c*d^3 + 49*d^4)*cos(f*x + e)^2 - (3*c^4 - 10*c^3*d - 24*c^2*d^2 - 54*c*d^3 + 85*d^4)*cos(f*x + e))
*sin(f*x + e))*sqrt(a*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c))/((a^3*c^5*d - 3*a^3*c^4*d^2 + 2*a^3*c^3*d^3
+ 2*a^3*c^2*d^4 - 3*a^3*c*d^5 + a^3*d^6)*f*cos(f*x + e)^4 - (a^3*c^6 - a^3*c^5*d - 4*a^3*c^4*d^2 + 6*a^3*c^3*d
^3 + a^3*c^2*d^4 - 5*a^3*c*d^5 + 2*a^3*d^6)*f*cos(f*x + e)^3 - (3*a^3*c^6 - 4*a^3*c^5*d - 9*a^3*c^4*d^2 + 16*a
^3*c^3*d^3 + a^3*c^2*d^4 - 12*a^3*c*d^5 + 5*a^3*d^6)*f*cos(f*x + e)^2 + 2*(a^3*c^6 - 2*a^3*c^5*d - a^3*c^4*d^2
 + 4*a^3*c^3*d^3 - a^3*c^2*d^4 - 2*a^3*c*d^5 + a^3*d^6)*f*cos(f*x + e) + 4*(a^3*c^6 - 2*a^3*c^5*d - a^3*c^4*d^
2 + 4*a^3*c^3*d^3 - a^3*c^2*d^4 - 2*a^3*c*d^5 + a^3*d^6)*f - ((a^3*c^5*d - 3*a^3*c^4*d^2 + 2*a^3*c^3*d^3 + 2*a
^3*c^2*d^4 - 3*a^3*c*d^5 + a^3*d^6)*f*cos(f*x + e)^3 + (a^3*c^6 - 7*a^3*c^4*d^2 + 8*a^3*c^3*d^3 + 3*a^3*c^2*d^
4 - 8*a^3*c*d^5 + 3*a^3*d^6)*f*cos(f*x + e)^2 - 2*(a^3*c^6 - 2*a^3*c^5*d - a^3*c^4*d^2 + 4*a^3*c^3*d^3 - a^3*c
^2*d^4 - 2*a^3*c*d^5 + a^3*d^6)*f*cos(f*x + e) - 4*(a^3*c^6 - 2*a^3*c^5*d - a^3*c^4*d^2 + 4*a^3*c^3*d^3 - a^3*
c^2*d^4 - 2*a^3*c*d^5 + a^3*d^6)*f)*sin(f*x + e))]

Sympy [F(-1)]

Timed out. \[ \int \frac {1}{(3+3 \sin (e+f x))^{5/2} (c+d \sin (e+f x))^{3/2}} \, dx=\text {Timed out} \]

[In]

integrate(1/(a+a*sin(f*x+e))**(5/2)/(c+d*sin(f*x+e))**(3/2),x)

[Out]

Timed out

Maxima [F]

\[ \int \frac {1}{(3+3 \sin (e+f x))^{5/2} (c+d \sin (e+f x))^{3/2}} \, dx=\int { \frac {1}{{\left (a \sin \left (f x + e\right ) + a\right )}^{\frac {5}{2}} {\left (d \sin \left (f x + e\right ) + c\right )}^{\frac {3}{2}}} \,d x } \]

[In]

integrate(1/(a+a*sin(f*x+e))^(5/2)/(c+d*sin(f*x+e))^(3/2),x, algorithm="maxima")

[Out]

integrate(1/((a*sin(f*x + e) + a)^(5/2)*(d*sin(f*x + e) + c)^(3/2)), x)

Giac [F]

\[ \int \frac {1}{(3+3 \sin (e+f x))^{5/2} (c+d \sin (e+f x))^{3/2}} \, dx=\int { \frac {1}{{\left (a \sin \left (f x + e\right ) + a\right )}^{\frac {5}{2}} {\left (d \sin \left (f x + e\right ) + c\right )}^{\frac {3}{2}}} \,d x } \]

[In]

integrate(1/(a+a*sin(f*x+e))^(5/2)/(c+d*sin(f*x+e))^(3/2),x, algorithm="giac")

[Out]

integrate(1/((a*sin(f*x + e) + a)^(5/2)*(d*sin(f*x + e) + c)^(3/2)), x)

Mupad [F(-1)]

Timed out. \[ \int \frac {1}{(3+3 \sin (e+f x))^{5/2} (c+d \sin (e+f x))^{3/2}} \, dx=\int \frac {1}{{\left (a+a\,\sin \left (e+f\,x\right )\right )}^{5/2}\,{\left (c+d\,\sin \left (e+f\,x\right )\right )}^{3/2}} \,d x \]

[In]

int(1/((a + a*sin(e + f*x))^(5/2)*(c + d*sin(e + f*x))^(3/2)),x)

[Out]

int(1/((a + a*sin(e + f*x))^(5/2)*(c + d*sin(e + f*x))^(3/2)), x)